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Abstract 

This is a review of the recent progress on a holographic description of the 
Schwinger effect. In 2011, Semenoff and Zarembo proposed a scenario to study the 
Schwinger effect in the context of the AdS/CFT correspondence. The production 
rate of quark anti-quark pairs was computed in the Coulomb phase. In particular, it 
provided the critical value of external electric field, above which particles are freely 
created and the vacuum decays catastrophically. Then the potential analysis in the 
holographic approach was invented and it enabled us to study the Schwinger effect 
in the confining phase as well. A remarkable feature of the Schwinger effect in the 
confining phase is to exhibit another kind of the critical value, below which the pair 
production cannot occur and the vacuum of the system is non-perturbatively stable. 
The critical value is tantamount to the confining string tension. We computed the 
pair production rate numerically and introduced new exponents associated with the 
critical electric fields. 
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1 Introduction 


The Schwinger effect is apair creation process in quantum electrodynamics (QED), 


due to an external electric hcl 
by 


The production rate (per unit time and volume) is given 


exp 


irm 


eE 


( 1 . 1 ) 


where m, e and E are an electron mass, an elementary electric charge and an external 
electric held, respectively. This formula is computed under the weak electric-held (eE <C 
m 2 ). A typical value of E for which the Schwinger effect becomes significant is estimated 
as 

( 1 . 2 ) 


m 


E = E Sch = — ~ 1.3 x 10 18 V/rn. 

e 


1 The Schwinger effect is not intrinsic to QED, but it is ubiquitous in quantum field theories coupled 
with a U{ 1) gauge field. Furthermore, it may be generalized to non-abelian gauge fields like color fields 
in quantum chromodynamics (QCD) [UjB]- 
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Thus -Esch is extremely large in comparison to typical values of electric field in table-top 
experiments. For example, a typical value of electric field necessary to ionize an atom is 
given by 

-^ion = m 2 e 5 Qig ~ 5.2 x 10 11 V/m, (1.3) 

with the fine-structure constant a s = e 2 /An. At least so far, in real experiments, the 
Schwinger effect has not been observed yet, but it may be observed in the near future. 
The XFEL project at DESY and the ELI project in Europe plan to produce extremely 
strong electric field very close to E Sch . 

The exponential factor in the production rate indicates that the pair production pro¬ 
cess should be described as a non-perturbative phenomenon like a tunneling process in 
quantum mechanics. To make a virtual pair of electron and positron be real particles, it 
is necessary to achieve larger energy than the static energy from an external source. A 
phenomenological potential, 

V tot (x) = 2m- — - eEx , (1.4) 

x 

leads us to a tunneling picture. The potential V tot (x) is composed of three parts: 1) the 
static electron mass, 2) the Coulomb potential with the distance x between electron and 
positron, and 3) an energy provided by the external electric field E. For the detail, see 
e.g. chapter 13 in the book [Tj. 

The potential shapes are depicted in Fig. |T|. When E is not so large, one can see the 
potential barrier (the blue line in Fig. [[]). But the particle can penetrate it quantum 
mechanically, the particle creation can be described as a tunneling phenomenon. This 
is nothing but the Schwinger effect. It would be worth estimating the production rate 
roughly in the WKB sense. The triangle approximation of the potential barrier leads to 
the following exponential factor, 



This is very close to Schwinger’s result ('ll.Ill , despite the rough estimation. 

As E increases, the barrier tends to be lowered and at last it disappears completely at 
a certain value of E (the red line in Fig.Q]). This value is referred to as the critical electric 
field E c . The disappearance of the barrier indicates that the QED vacuum becomes 
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Vtot{x) 



Figure 1: The potential shapes. For the blue line, the particle production is described as 
a tunneling process. The red line indicates that the potential barrier vanishes at a certain 
value E = E c (critical electric field). For the green line, the vacuum is catastrophically 
unstable. 

unstable catastrophically once E has reached E c . It is easy to evaluate this value as 

9 

„ rn _ 9 

eE c = — ~ 137m 2 . 
a s 

Note that this value is far beyond the weak-field approximation eE <C m 2 . Hence, frankly 
speaking, it has not been clarified whether the catastrophic decay would really occur or 
not. 

An intriguing question to be answered is “Does the catastrophic vacuum decay in 
QED really occur?” Although we have discussed the Schwinger effect in QED so far, let 
us change the direction a little here and focus upon a critical behavior in string theory. 
It is well known that there exists an upper critical value of electric field in the context of 
string theory m- This bound obeys from a self-consistency condition of string theory. 
On the other hand, we know intimate relations between string theories and gauge theories. 
For example, these are established through the AdS/CFT duality [T0l - fT2] . According to 
this duality, a string theory on an anti-de Sitter (AdS) space is equivalent to a conformal 
field theory (CFT). The most well-studied example is the duality between type IIB string 
theory on AdSs x S 5 and the AT = 4 supersymmetric Yang-Mills theory (SYM) with an 
SU(N ) gauge group in four dimensions. 

Thus, one may expect a connection of electric fields between the string theory and 
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gauge theory sides through the AdS/CFT duality. Then the result in the string-theory 
side suggests that the catastrophic vacuum decay can really occur at least in a class 
of gauge theories which have gravity duals. Thus, to investigate the possibility of the 
catastrophic vacuum decay in gauge theories much deeper beyond the weak-held approx¬ 
imation, it would be worth studying the Schwinger effect in the context of the AdS/CFT 
correspondence ra. 

Another motivation is to argue the Schwinger effect in QCD. It may give rise to a 
new mechanism of a confinement/deconfinement phase transition and it would have a 
connection with the RHIC and LHC experiments, where strong electro-magnetic fields 
and color fields are induced due to the collision of heavy ions. It motivates us to study 
the Schwinger effect in confining gauge theories. However, it is not easy to tackle this issue 
with the standard method in quantum field theories. A nice way is to employ a holographic 
computation by realizing confining gauge theories with appropriate D-brane set-ups. In 
fact, we have considered the Schwinger effect in confining gauge theories along this line 
[T4UT6] . The potential analysis [14f[l5j has been done for general confining backgrounds 
[IMH] by generalizing the procedure in the Coulomb phase ra- The production rate has 
been evaluated numerically in the case of a confining D3-brane background [16] (which is 
an example of AdS solitons M). 

The organization of this review is as follows. Section 2 gives a short review of the 
world-line instanton method based on articles j22l - l24] . This is a method to calculate the 
production rate and becomes a key ingredient in preparing a holographic computation 
in the subsequent sections. In section 3, we consider the Schwinger effect in the J\f — 4 
SYM. This system is conformal and contains no fundamental matter held. Hence, naively, 
the Schwinger effect cannot be argued. A possible way of overcoming these points is 
to employ the Higgs mechanism. We first compute the production rate in the Higgsed 
J\f — 4 SYM naively by using the world-line instanton method. After all, this computation 
leads to a puzzle on the critical value of external electric held. Then, to resolve this 
puzzle, an improved holographic set-up is presented by following the seminal work by 
Semenoff and Zarembo The production rate is reconsidered along this line and 

the puzzle is certainly resolved. Section 4 is the main part of this article based on a 
series of our works |[T4HT6] and investigates the Schwinger effect in a confining gauge 
theory by employing a confining D3-brane background [2T] . A remarkable feature is 
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that there is another kind of critical value of electric field due to the existence of the 


confining phase. We have evaluated the production rate numerically and computed new 
exponents associated with the critical electric fields. Section 5 is devoted to conclusion 
and discussion. 

2 The world-line instanton method 

In this section, we will give a short review of the world-line instanton method by following 
the works fZ2 \fM\ . 

2.1 The production rate at weak coupling 

There are various methods to calculate the pair production rate of the Schwinger effect 
in QED. We shall introduce here one of them, called the world-line instanton method. 
An advantage of this method is that one can easily add a contribution coming from the 
Coulomb potential. In particular, it enables us to compute the production rate at arbitrary 
coupling. In the following, we will explain first the computation at weak coupling, then 
the one at arbitrary coupling. 

As a simple example, let us first consider a massive scalar QED. The action for a 
massive scalar QED (in the Euclidean signature) is 



( 2 . 1 ) 


with the covariant derivative, 


( 2 . 2 ) 


= d lt + ieA^ + ieA™ . 


Note here that the vector field is divided into a dynamical (fluctuating) part A M and an 
external part A e *. The path integration is performed only for the dynamical part, not 
for the external part. 

The production rate T is given by the imaginary part of the vacuum energy density 


£ o, 



(2.3) 
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Here V 4 is the volume of the four dimensional space. By performing the path integral for 
the scalar held 0, the following expression is obtained, 


V 4 r = —2 Im log / VA exp 




tr log (— D 2 + m 2 )^ . (2. 


4 ) 


In a small coupling region, the dynamical part of the vector held can be ignored. Hence 
the expression (12.41) is approximated as 


VaT ~ —2 Im 


'0 


’d T 
~T 


Vx exp 


2T 


dr x~ 


m 2 T 


+ ie dr A^x 


’o 


n m 


(2.5) 


Here we have used Schwinger’s parametrization 

r dt 

log a = - / — e 

Jo 1 


—at 


and the quantum-mechanical path-integral representation. Note that x^{r) satishes the 
periodic boundary condition r M (0) = x M (l). 

Next, as for the expression (12. 5 j) , the T-integration is hrstly performed, and then the 
integral for x(r) is donc^l. For later convenience, let us suppose that 


m\ 


dr x 2 1. 


( 2 . 6 ) 


This condition (12.6(1 will be identified with the weak-held condition for the external electric 
held, as we will see later. It is worth noting that the integration about T can be regarded 
as a modified Bessel function, 


dt 


K o(x) = -l —exp (-t- — 


t 


x 


and its asymptotic behavior for large x is 


K 0 (x) ~ 





2r 


Then, under the condition (12.6(1 . the expression ((2.5(1 is simplified as 


V 4 T = — 2 Im / Vx ■ 
where particle an d T 0 are dehned as 

^particle = Tfl 


2vr 



mV Tr 


exp ( Fpartide) , 


(2.7) 


( 2 . 8 ) 


dr x 1 — ie / dr x^ , 


( 2 . 9 ) 


The integrations in the opposite order lead to the Heisenberg-Eular Lagrangian [2]. 
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( 2 . 10 ) 


To = 



The above computation is identical with the method of steepest descent for the T- 
integration. But, note that the condition (12.61) is not necessary if we do not want to 
use the asymptotic form of the Bessel function at this stage. 

The path integral about x{r) is evaluated by the method of steepest descent. The 
equation of motion obtained from S^articie is 


mXfj, 

yjl dT±2 


—ieF^x 


V • 


Multiplying (I2.11|) by and integrating it, we find the following relation, 


( 2 . 11 ) 


x 2 = a 2 (constant). 


( 2 . 12 ) 


Let us suppose a constant electric field is turned on the aq-direction. Then the vector 
potential is taken as 


Af^o) = ~ iEx o ■ 


(2.13) 


The other components are set to be zero. Note that the imaginary unit i appears because 
we are now working in the Euclidean signature. 

With the constant electric field, the classical solutions are represented by 


m 

x 0 = — cos 
eE 


( aeE 

V m 


-T 


m 


x\ = —— sm 
eE 


. (aeE 


\ m 


-T 


X2 = X 3 = 0 . 


(2.14) 


The periodic boundary condition x M (0) = x^(l) determines a like 

aeE 


= 2nn , 


n e N. 


(2.15) 


m 

Thus the classical motion has been completely fixed. This solution is often called “in- 
stanton,” just because this is a classical solution in the Euclidean signature. 

With the relation (12.151) . it is easy to check that the assumption (12.61) is nothing but 
the weak-field condition 


E E Sch , 


(2.16) 


as noted before. Furthermore, the assumption is the same as the condition to ensure that 
the steepest-descent method is a good approximation. 
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By combining the classical action 

q(n) 

Wl 


7i m 2 n 
eE 


with one-loop prefactors, the production rate is evaluated as 


r = 


[eEf 
(2vr) 3 


E 

n=l 


-i) 


n +1 


n z 


■ exp 


irm 


eE 


-n 


(2.17) 


The production rate (I2.17P is the same as the Schwinger formula for the scalar QED. This 
expression is valid at small coupling and under the weak-held condition. The method 
presented here is applicable to a spinor QED as well [23]. 

We have not explained the one-loop prefactor here. For the detail of the derivation, 
for example, see the works [ 221 [24] . 


2.2 The production rate at arbitrary coupling 


Let us generalize the production rate (12.17ft at weak coupling to the one at arbitrary 
coupling, though we still suppose the weak-held condition. 

Before going to the detail, it is worth to see a heuristic argument based on the pair 
production rate of monopole and anti-monopole. In the Georgi-Glashow model, it is 
evaluated as [25] 

(nR \ 2 / ttM 2 n 2 \ 

(2.18) 


r (9-S) 2 

r = exp 


(2tt) 


t rM 2 g 2 

-b — 

gB 4 


where B, M = Airm/e 2 and g = 47r/e are an external magnetic held, a monopolc mass 
and a magnetic charge, respectively. This result (12.18ft is valid under the following two 
conditions: 

g 2 > 1, gB <C M 2 . 

Let us perform electric-magnetic duality for the above result. The production rate is 
mapped to 

(f 'Trm^ 

(2.19) 


„ (e-E) 2 

T - — .. exp 


(2x) 


2 2 
7rm e 

~^E + T 


But the validity region is 


e 2 > 1, eE < m 2 , 


and the strong coupling region is supposed as opposed to the previous computation. This 
observation indicates the expression (I2.19P should be valid for arbitrary coupling. In fact, 
this expression can be obtained by a direct computation [22], as we will show below. 


















In the previous subsection, the coupling constant has been supposed to be so small 
that the dynamical held A M is ignorable. But we will take account of at finite coupling 
here. Then the equation (j2.4j) is replaced by 

V 4 T = —2 Im log (exp [—tr log (— (d + ieA + ieA ex ) 2 + m 2 )]) , (2.20) 


where the expectation value is defined as 


(9(A)) 



Furthermore, the equation (12.201) can be approximated as 


U 4 r ~ 2 Im (tr log (— (d + ieA + ieA ex ) 2 + m 2 )) 


( 2 . 21 ) 


( 2 . 22 ) 


under the weak-held condition. This can be regarded as an electron loop expansion in the 
Feynman diagrams. 


Here we keep the diagrams with a single electron loop and drop off the higher-loop 
contributions. This selection is a good approximation when the external held is weak 
enough. By doing the same calculation as in the previous subsection, the production rate 
is obtained as 


f°° d T f 

V 4 T = —2 Im / — / Dx exp 


x 


T 


exp \ie ® An da; 


-/ d t x 2 

2TJ 0 


m 2 T 


+ ie(j) A ex dx A 


(2.23) 


In comparison to (12.51) . the expression (12.231) includes the contribution coming from the 
Coulomb interaction. This modification is rehected as a 17(1) Wilson loop in (I2.23F 
Under the condition (12.61) , the integration about T leads to 

W * -2 In, jv* exp (-, W - £ // ) , (2-24) 

where we have employed the following identity: 


exp i le 


A, 




dx ■ dy \ 
(x-y) 2 ) ' 


(2.25) 


The derivation of (12.251) is straightforward. 
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Let us next perform the path-integral about x(r ). The Coulomb interaction should 
change classical solutions, as a matter of course. However, note here that the Coulomb 
interaction term is invariant under a scale transformation and a rotation in the plane on 
which the instanton is rotating. Hence it does not depend on the size of the classical 
solutions. Thus the classical solutions are not changed when the Wilson loop is added as 
a perturbation. 

The classical action is modified as 


Q (n) _ 
°cl — 



Here we have used the formula 



dx ■ d y e 2 

7-WT = n - 7 

(x — y) 2 4 


(2.26) 


(2.27) 


where an unphysical divergence has been ignored. 

Recall that the contributions of the Coulomb interaction are small. Then the prefactor 
of the exponential has not been modified. As a result, the production rate is given by 


T 


( e £) 2 y 

( 2 7t) 3 


(-l ) n+1 

T) 2 


exp 



(2.28) 


Thus we have reproduced the conjectured form (12.191) as the n — 1 case. 

Note that the production rate is not exponentially suppressed any more just after E 
has reached the critical value 


E = Er = 


4:7im 2 


That is, the vacuum becomes unstable catastrophically above this critical value. However, 
the weak-held condition E <C i?sch is broken as follows: 


£c = — =► E c ~137E Sch . 

oi s 

Hence it is not certain whether the catastrophic vacuum decay really occurs or not. ft is 
significant to answer this question. A possible way is to follow a holographic computation, 
as we will introduce in the next section. 
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3 The Schwinger effect in the AdS/CFT correspon¬ 
dence 

In this section, we will consider the Schwinger effect in the context of the AdS/CFT 
correspondence. There are many variations of AdS/CFT, but we will concentrate on the 
most typical example, the duality between type IIB superstring on AdS 5 x S 3 * 5 and the 
J\f — 4 SU(N ) SYM in four dimensions. However, in order to argue the Schwinger effect, 
there are three obstacles: 

1. The A f = 4 SU(N) SYM is conformal (i.e., all of the fields are massless). 

2. All of the matter fields belong to the adjoint representation. 

3. There is no 17(1) gauge field. 

As for (3), as a matter of course, one may consider non-abelian Schwinger effects. But 
it is not so easy to study the issue itself. Hence it is better to focus upon an abelian 
Schwinger effect as a first trial. 

Let us first introduce a 17(1) gauge held and fundamental matter fields by employing 
the Higgs mechanism. Then one can compute the pair production rate of the fundamental 
matter fields by following the procedure introduced in Sec. 2. However, we will encounter a 
puzzle for the critical value of electric held. To resolve it, Semenoff and Zarembo proposed 
a prescription and presented an improved holographic set-uj . Subsection 3.2 will 
introduce their proposal. 

3.1 Set-up 

We start from the J\f —4 SU(N + 1) SYM theory. It consists of a gauge held A M (p = 
0, • • • , 3), six real scalar helds (/ = 1, ■ ■ ■ ,6) and four Weyl fermions T, the hat is 
attached for later convenience. All of the helds belong to an adjoint representation of 
SU(N + 1). One should notice here that the A f= 4 SYM includes neither a 1/(1) gauge 
held and nor fundamental matters. 

3 For an earlier trial, see the work [26]. In Ref. m , it has been shown that an electric field creates a 

pair production current on a probe brane at zero temperature. The finite-temperature case is discussed 

in Refs. [MIES]. 
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Then, by breaking SU(N + 1) to SU(N) x U( 1) with the Higgs mechanism, let us 
introduce a 17(1) gauge field. The SU(N + 1) fields are decomposed into the SU(N) part, 
the U( 1) part and non-diagonal parts: 


A 







(3.1) 


Here A M [a^], [ 0 /], and T \ip\ are the SU(N) [t/(1)] gauge, the scalar and the fermionic 
fields, respectively. The vacuum expectation value (VEV) of the scalar fields is supposed 
to be 

6 

($/) = diag(0, ■ ■ ■ , 0, m0/), 5^0/= 1. 

7=1 

In terms of D-branes, the Higgs mechanism corresponds to separating a single D3-brane 
from the remaining stuck of N parallel D3-branes. The distance between the separated 
D3-brane and the N D3-branes is related to the VEV of the scalar fields. 

As a result, the J\f —4 SU(N + 1) SYM action *S^^ym i s decomposed like 

q SU(N+ 1) _ qSU(N) q U( 1) Q 

'VV’=4 SYM ~ *AV=4SYM + Aa/=4SYM + ‘-’W • 


Here Sw is the action for the W-boson supermultiplet 
1 


S'w — 


9ym 


d x 


D^ujj + l o\ (<E>x — md K ) 2 ojj — m 2 u\9i9jUj + • • • 


(3.2) 


For our purpose, we will concentrate on the complex scalar fields u>i (called “quarks”) 
and drop off the vector field and the fermionic field x i n the W-boson supermultiplet. 
The higher order interactions are also ignored. 

A remarkable point is that the covariant derivative D fl in (13.2ft is given by 


Df,. = d fl + ia M - iA t 


Thus the scalar fields in the W-boson supermultiplet have a finite mass m and couple to 
the U( 1) gauge field as well as the SU(N) gauge field A M . 

Note here that five scalar fields are massive but one scalar field is massless. For 
example, let us consider the case with 9j = (0, 0, 0, 0, 0,1). For I — 1,... ,5, mass terms 
exist in the action (13.2ft , but for / = 6 it vanishes. The massless field is absorbed into 
the vector field u as the longitudinal mode so that becomes massive. 
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Let us calculate the pair production ratqj for ujj by using the world-line instanton 
method. In the following, we take the ’t Hooft limit [30] where N —> oo with A = Dym-^ 
fixed. The value of A is also supposed to be sufficiently large so as to suppress the 
dynamical part of the 1/(1) gauge field. Then can be regarded as an external field like 


a® x = —iEx o 


(the other components are zero). 


(3.3) 


Note that we work in the Euclidean signature to calculate the production rate. 

Then the production rate can be evaluated aj^| 

E 4 T = — 5Nlm jvxg[x{T)] exp m J d J d Ta e *x^\ (W[x \), (3.4) 

dr (iA^x^ +\ , (3.5) 

where g[x(r)\ is a functional of x(t) and does not contribute to the exponential factor in 
the production rate. In the derivation of (13.41) . we have assumed that quark mass is very 
heavy, and the external electric field is weak, E -C m 2 . 

It would be helpful to comment on some differences between (12. 8ft and (13.41) . Since 
the production rate is proportional to the number of fundamental particles, 5 N appears 
in (13.4[) (note that one of loj has been eaten by oj^). The kinetic term for the instanton 
is changed, but this change does not affect the exponential factor in the leading-order 
computation. Note also that W[x] in (13. 4 j) is a supersymmetric SU(N) Wilson loop [341135] 
because the SU(N) gauge field A M has been regarded as a dynamical field. 

Now it is necessary to evaluate the VEV of the Wilson loop. One can compute it by 
examining the area of the minimal surface of a string attaching to the circle. The result 
at strong coupling is given by [ 331(36] 


W[x] = tr su(N) exp 


(W[x\) =e^\ 

4 Recall that w/ , and \ have the same mass. Hence u>i , and x may be produced equally. In 
fact, the exponential factor does not depend on spins, though the prefactors are different. In addition, as 
for , there is another kind of problem called the Nielsen-Olsen instability. When pairs are created 
under external electric and magnetic fields. See the work [31: for the detail of this instability in the 
holographic Schwinger effect. 

5 For the detailed derivation, see sec. 2 in the work [32]. It is basically the same as the derivation of a 
supersymmetric Wilson loop in an earlier work [33], up to prefactors. 
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by supposing that the mass of quarks should be very heavy. Then the exponential factor 
in the production rate is evaluated as 

+ y/\') . (3.6) 


T ~ e 5cl = 


exp 


Tim 

~E 


This is an expected form because \f\ is regarded as a coupling constant in the large N 
gauge theory |[34jj35]. From (13.6ft . one can read off the critical electric field, 

.2 


7T m 


(3.7) 


When E < E c , the pair production is suppressed exponentially. When E > E c , the 
vacuum becomes unstable catastrophically. 

ft may seem that the critical electric field (13.711 satisfies the weak-field condition be¬ 
cause A is supposed to be very large. However, it would be unlikely that heavy particles 
are produced due to the Schwinger effect because such a process is severely suppressed. 
Hence the above computation may be unsatisfactory. 


In addition, note that the critical electric field (13. 7 p does not agree with the one 
obtained from the DB1 action of a D3-brane sitting near the AdS boundary^, 

2n m 2 


E DBI = 




(3.8) 


This disagreement may lead to a puzzle. 

In the next subsection, we will introduce a nice prescription to resolve the above two 
points simultaneously. 


3.2 Semenoff-Zarembo’s prescription 

In the work [13], Semenoff and Zarembo proposed a holographic set-up to study the 
Schwinger effect in the higgsed J\f — 4 SYM. 

The set-up is the following. A single D3-brane is separated from the stuck of N 
parallel D3-branes. The near-horizon limit of the stuck of N D3-branes is represented by 
the AdSs x S 5 geometry: 

ds 2 = gMN<ix M dx N 

6 If the D3-brane is located at the boundary, then the mass m diverges. Hence m has become finite 
by slightly separating the D3-brane from the boundary. 
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Boundary 


00 



Probe D3-brane 



r o 


-L r = 0 

Horizon 


Figure 2: The location of the probe D3-brane and the configuration of the string world-sheet. 

= ^—dx^dx^ + ^5-dr 2 + L 2 dfi 2 , (3.9) 

L z r z 

where L = A is the curvature radius. The coordinates x M (M = 0,..., 9) describe 

the ten-dimensional spacetime and x^ (/j, = 0,... ,3) represent the four-dimensional space- 
time in which the dual gauge theory lives. The coordinate r is the radial direction of the 
AdS space. The horizon of AdS is located at r = 0 and the conformal boundary is at 
r = 00 . The separated D3-brane can be treated as a probe brane in the bulk AdS. A 
remarkable point is that the probe D3-brane is put at an intermediate position (r = ro) 
between the horizon and the boundary, and it sits at a point on S' 5 . For the configuration, 
see Fig. [2j 

According to the set-up, the mass of a single quark is not infinite any more but finite. 
The quark mass is measured by the energy of a string stretched from the probe D3-brane 
to the horizon like 

m = T F r 0 . (3.10) 

Here Tp = 1/27T a' is the string tension. 

Inspiring from the fact that a circular Wilson loop appears in the production rate 
113. 4 jl . Semenoff and Zarembo proposed a new scenario that the production rate may be 
computed by a circular Wilson loop on the probe D3-brane equipped with a constant 
electric H-held. The VEV of the Wilson loop can be evaluated by the area of a minimal 
surface of the fundamental string attaching to the boundary loop. That is, the exponential 
factor in the production rate is estimated as 

T ~ exp(—SW g - S Ba ) ■ (3.11) 
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Here Sng and S Bs are the Nambu-Goto (NG) action and the coupling to an NS-NS 2-form. 
In the Euclidean signature, these quantities are given by, respectively, 


Sng — T f J d 2 cr y /det G a p , (3-12) 

S Bs = -T F J d 2 aB M Nd T x M d a x N . (3.13) 

The string world-sheet is parametrized by a a = (r, a). Then G Qj g is the induced metric 
and Bmn is an anti-symmetric two-form flux. Note again that we work in the Euclidean 
signature for both the bulk spacetime and the world sheet. 

It would be worth comparing the world-line instanton method and Semenoff-Zarembo’s 
prescription. Recall that in the former method the following expression 

m J dr — log(W[x]) + i Jdra^x^ (3-14) 

appears in the exponential in (I3.4[) . In the latter method, the factor (13.14}) is replaced by 
Sng + Sb s ■ That is, the particle action and the VEV of the Wilson loop are replaced by 
the NG string action Sng- Then the NS-NS two-form in S B2 is interpreted as an external 
electric field on the probe D3-brane via the following relation, 

TfBqi = E. 


In the ’t Hooft limit, the classical analysis is sufficient on the string-theory side. To 
begin with, we construct a classical solution whose boundary is a circular Wilson loop on 
the probe D3-brane by taking the following ansatz (See Fig. [2]) : 


x = x(a) cos t , 




r = r(a ), 


x 


x = x(a) sm r , 
2 \ 2 

= R 2 -f 


r 0 

icfl. 


x 2 = x A = 0 , 


, (3.15) 

The range of r is 0 < r < 27r. 


where R is the radius of the circle on the probe D3-brane 
Note that the coupling to a constant NS-NS 2-form does not change the classical equation 
of motion. Hence the classical solution is given by 


X(<7 = 


cosher 


'R 2 + I —) 
To J 


Her = 


tanh a 


(3.16) 


7 Note that a; 0 = x (a) cos (nr), x 1 = x(cr) sin(nr) with n > 2 satisfy the circular ansatz (111.151) . 
However, the contributions of the solutions with n > 2 are suppressed in comparison to the one with 


n = 1, and hence we will concentrate on the ansatz (13.151) . 
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Here the range of a is taken as <7o < cr < oo and sink cr 0 = E 2 / Rr$ . 

The next is to argue the boundary condition on the probe D3-brane. Note that the 
NS-NS 2-form is relevant to the boundary condition, while it is blind to the equations of 
motion. In the present case, one needs to impose the following mixed boundary condition: 

f) r M a M 

VfetGG^gMNj— - B mn — = 0 . (3.17) 

oa a or 


This boundary condition determines the radius R like 


R = 



(3.18) 


Here we have defined E c as 




2nm 2 


(3.19) 


Note that the classical solution does not exist when E is larger than E c . It is also possible 
to derive the relation (13.18[) by taking a variation of S c \ with respect to R , as originally 
argued in the work [15]. 

Finally, by substituting the classical solution (I3.16P to the string action, the production 
rate is evaluated as 


T ~ exp 



= exp 


nm 2 

~JT 


+ Vx 



(3.20) 


The first and second terms in the exponential are the same as the expression (13.61) . The 
third term should be regarded as a correction term and can be ignored under the weak-held 
condition E -C m 2 . In fact, the classical action vanishes when E = E c . 

It is worth comparing the result (13.19j) with the critical electric field derived from the 
DBI action of the probe D3-brane. Here we will work in the Lorentzian signature. The 
DBI action in the bulk AdS is given by 


SnRT — — T] 


D3 


L 4 


(TX \ 1 — 


(2na') 2 L i 


E 2 


where the D3-brane tension is given by 

2d3 = 


g s (2n) 3 a 


12 


(3.21) 
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One can see that the DBI action is ill-defined when E > 2irm 2 /\/X. Thus, the critical 
electric field is given by 

£ c dbi = ^ . (3.22) 

This completely agrees with the result (13.19j) . 

As we have seen so far, Semenoff and Zarembo’s proposal works well on the critical 
electric field and has resolved the two unsatisfactory points raised at the end of sec. 3.1. 
Then the production rate (I3.20j) coincides with the result (13.61) computed with the world¬ 
line instanton method, under the weak-held condition E <C m 2 . Thus it seems likely that 
there would be no problem for the exponential factor. 

However, it has not succeeded yet to evaluate the prefactor of the exponential with 
the present holographic set-up. The difficulty comes from the fact there is a subtlety 
in regularizing quantum fluctuations around a circular Wilson loop, in comparison to a 
straight line. In fact, some trials have already been done [371 - HO] . but unfortunately none 
of them has succeeded. 

Finally, it is useful to comment on some generalizations of Semenoff and Zarembo’s 
work [13J. One may consider the finite temperature case |20lf3lj by following the works 
gum]. Bolognesi et. al. ED argued the pair creation of monopolc and anti-monopole 
in an external magnetic field. Some cases with external electro-magnetic fields are also 
studied [32]. 

4 The Schwinger effect in a confining gauge theory 

So far, we have studied a holographic description of the Schwinger effect in the Coulomb 
phase. An interesting direction is to consider the Schwinger effect in the confining phase 
by generalizing the previous set-up. This section is devoted to the study of the Schwinger 
effect in a confining gauge theory by employing a holographic approach mm- 

It would be important to motivate the readers to study the Schwinger effect in the 
confining phase, before going to the detail. Recently, QCD in strong external fields 
has been intensively studied to explain experimental data of RHIC and LHC. In these 
experiments, a extremely strong magnetic field is induced, because colliding nuclei have 
large electric charge and move very fast. This strong magnetic field may have a visible 
effect on the nuclear dynamics. Apart from the magnetic field, a strong electric field is also 
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created around the nuclei as well as a strong color field. Thus, to capture the underlying 
physics concerned with the experiments, it is inevitable to study the QCD dynamics in 
the presence of strong external fields. As a matter of course, the Schwinger effect in the 
confining phase is included in this direction. 

We will utilize a holographic approach to study the Schwinger effect, but one may 
wonder whether the lattice formulation would be enough to study it. In the presence of 
an external electric field, one encounter a notorious sign problem and hence it is quite 
difficult to use the lattice formulation in general. Some researchers may hesitate to bor¬ 
row a holographic computation. But one cannot disguise the fact that there is no definite 
method to tackle this issue instead of it, at least so far. It would be quite useful to 
capture qualitative pictures of the Schwinger effect in the confining phase even by em¬ 
ploying the AdS/CFT. In fact, one can see quantitative predictions as well as qualitative 
understanding by pursuing this direction. 


4.1 Set-up 


First of all, let us introduce the set-up to study the Schwinger effect in the confining phase. 
For this purpose, the bulk AdSs has to be replaced with another confining background 
with a dimensionful parameter. There are lots of backgrounds dual to confining gauge 
theories. A simple example is an AdSs soliton background (often called a confining D3- 
brane background)^. 

The metric of the AdSs soliton background is given by 

2 


ds 2 = =- -(dx 0 ) 2 + V(d^) 2 + f(z)(dx 3 ) 2 + 
z z 

/(*) = !- /W> 


2 — 1 


dz 2 

W) 


+ L 2 dn 


5 5 


(4.1) 


Now the radial direction is described by z and a scalar function f(z) contains a constant 
parameter z t . Then the geometry is cut off at z = z t and the conformal boundary is at 
z = 0. The parameter zt has the dimension of length and the inverse has the dimension 
of energy. Hence l/z t gives rise to the confining string tension on the dual gauge-theory 
side. As z t —3 > oo, f(z) —> 1 and the usual AdS 5 geometry is reproduced. 

8 In this review, we will concentrate on an AdSs soliton background for simplicity, but it would be 
easy to generalize the present analysis to the other confining backgrounds by following the work nsj. 
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As another point, the x 3 -direction is compactihed on a circle S 1 with the radius R = 
TTZ t . Hence the limit z t —* oo corresponds to the decompatification limit R —> oo. 
Anyway, the dual gauge theory lives on M 1,2 x S 1 rather than M 1 ’ 3 . Thus the dual gauge 
theory cannot be regarded as the usual confining gauge theory in 1 + 3 dimensions. Still, 
however, it would be possible to extract qualitative behaviors by examining the AdS 5 
soliton background. The results obtained here would be a key ingredient towards the 
study in the real QCD. 

4.2 Potential analysis 

Let us first study the Schwinger effect in the confining phase from the viewpoint of poten¬ 
tial analysis. Naively, the Schwinger effect could be argued with the modified potential 

V tot (x) = 2m + V(x) - Ex + a st x . (4.2) 

Here the modification is that a confining potential with the string tension cr st is added as 
well as the usual potential analysis. An important observation is that the string tension 
cr st and the electric field E compete with each other. When <r s t > E , the potential diverges 
as x —> oo and hence the Schwinger effect would not occur. On the other hand, when 
cr st < E , the situation is not so different from the Coulomb phase and the qualitative 
behavior would be almost the same as the previous. 

In the holographic set-up, it is necessary to evaluate the VEV of a rectangular Wilson 
loop for the quark anti-quark potential, ft can be computed with the minimal surface 
of a string attaching to the loop on the probe D3-brane embedded in the AdSs soliton 
background. 

In this subsection, we will work with 

_ L 2 
z 

as the radial coordinate, following the notation of our article p3]. Hereafter, we move to 
the Euclidean signature. The probe D3-brane is supposed to be at r = r 0 . We impose 
the following ansatz for the rectangular Wilson loop (See Fig. |31): 

x° = r , x 1 = a, r = r(a). (4.3) 
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Figure 3: The configuration of the string world-sheet for the quark and anti-quark potential. 


The NG part of the string action is 

C = V det G = 


1 / d r\ r 4 


1 _ r ^/ r i ydcr/ L 4 

Since the Lagrangian does not depend on cr explicitly, the Hamiltonian 

dC 


d(d a r 


-d a r — C 


(4.4) 


(4.5) 


is conserved. When we impose the boundary condition at the tip of the minimal surface 
like 


dr 

da 


= 0 , r = r c (r t < r c < r 0 ), 


the conserved quantity becomes 


= const. = r 2 n L 2 


dr \ r 4 


(4.6) 


(4.7) 


1 — r 4 /r 4 \da / L 4 
By rewriting the conserved quantity, we obtain 




(4.8) 


By integrating the expression (14.8p . the distance between a quark and an anti-quark, x, 
is given by 

2 L 2 f l/a d y 


x = 


r o° J i \/{y A - 1 )(y 4 - ( b/a ) 4 ) ’ 


(4.9) 
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where we have defined dimensionless quantities as 


r r c r t 

y = —, a = — , b = — . 
r c r o r 0 

The sum of the potential energy and static energy is evaluated as 

,4 


y 


rx/2 1/a. 

^pe+se = 2T f / dcr jC = 2Tpro& / d?/ — . - 

to A 


(4.10) 


Here x is a function of a , and hence the potential is a function of x thorough a. 

For large x limit (i.e. a —> b limit), the sum of the potential energy and static energy 
behaves as 


bpE+SE — Tp ( b 2 x + 2Tpr 0 b I - — 1 


(4.11) 


The first term is the quark and anti-quark potential with a confining string tension 


= T F 


n 


and the second term is static mass of quark and anti-quark 

2T F (r 0 — r t ) = 2 m w . 

By including the energy coming from the external electric held, the total potential is 
rewritten into an integral representation, 


ktot — kpE+SE — Ex 

rl/a 

= 2T F r 0 a / d y — ---- 

Ji /(y 4 - 1 ){y A - (6/a) 4 ) 


y 


2T F r 0 a ^ a 


d y 


a J i /(|/ 4 - 1 ){y A - ( b/a ) 4 ) ’ 

where we have introduced a dimensionless electric held a 


a = 




(4.12) 


(4.13) 


E c ' “ 1 L 2 

E c is the critical electric held obtained from the DB1 action. In the following, we define 
a dimensionless electric held a normalized by E c . 

The total potential is plotted in Fig. [4] as a function of x . Figure 0] shows the existence 
of two critical electric helds. The hrst one is 


E = E s = cr st (a = b 2 ), 


(4.14) 
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Figure 4: The plots of the total potential with b = 0.5 and 2L 2 /ro = 2Tftq = 1. The blue line is for 
a = 0.1. There is no zero other than the origin and hence the Schwinger effect does not occur. The red 
line is for a = 0.25. This value corresponds to b 2 = 0.25 and the potential becomes flat as x —> oo. For 
the values of a between 0.25 and 1.0, the potential barrier is formed and the Schwinger effect can occur 
as a tunneling process. When a = 1.0 (the orange line), the barrier just vanishes and the system becomes 
unstable catastrophically. 


and the second is 

E = E c (a = 1). (4.15) 

When the electric held is smaller than the confining string tension, the pair production 
is prohibited because the potential does not dump at infinity. When the electric held is 
larger than it, it is allowed as a tunneling process. Thus, at the hrst critical electric held, 
a confinement/deconhnement transition starts. As a side remark, a similar behavior of 
the phenomenological potential in a confining theory has been argued by using a lattice 
formulation [43]. 

The second critical electric held is the same as the critical electric held obtained in 
the Coulomb phase. That is, the vacuum becomes unstable catastrophically above this 
value of the electric held. Note that this qualitative behavior can also be understood 
analytically. For the detail, see the works mmm- 

fiere we are mainly concerned with the long-distance behavior of the potential. But 
it would be worth commenting on the short-distance behavior. As shown in Fig.Q] the 
potential becomes zero at the origin, in comparison to the usual potential analysis. This 
is just because the Coulomb potential is modified on the probe D3-brane sitting at an 
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intermediate position in the bulk. The modified Coulomb potential was originally argued 
by Kabat and Lifschytz |44]. 


4.3 The production rate 

In this subsection, we will compute the production rate numerically and examine two 
critical electric fields. Both of the resulting critical values agree with the ones obtained 
in the potential analysis D31- Furthermore, we will introduce new exponents associated 
with the critical behaviors and evaluate the numerical values. 

To compute the production rate, one has to evaluate the VEV of a circular Wilson 
loop on the probe D3-brane. Then the D3-brane is placed at an intermediate position 
between z = z t and 0 = 0. 

The first is to construct a classical string solution attaching to the circular Wilson 
loop. Let us suppose the following ansatz: 

x° = x(a) cost , x l = x(cr) sinr, z = z(a). (4.16) 

The other components are set to be zero. Suppose that the world-sheet coordinates (r, a) 
are restricted to the following range: 

0 < r < 2tt , 0 <cr<cr 0 - 

Then boundary conditions for x(a) and z(c r) are imposed like 

x(0)=0, x(a 0 ) = R, z(0) = z c , z(a 0 ) = z 0 - (4.17) 

The configuration is depicted in Fig. [5] 

In addition, a constant field B 2 is set to be 


T-fBqi = E. 


This B 2 field induces an external electric field on the probe D3-brane. Then the resulting 
NG action and the coupling to B 2 are given by, respectively, 

r R 


S'ng — 27tL Tp 


np I r /2 

w 


Jo z- 

r R 

Sb 2 — —‘J'kTyBqi I dxx = —ttER 2 . 


(4.18) 

(4.19) 
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Figure 5: The configuration of the string world-sheet for the Schwinger pair production. 


Here we take that a = x(a) with a diffeomorphism invariance. The prime denotes the 
derivative with respect to x . 

. Then the equation of motion for z(x) is obtained as 


, 2 xf(z) „ 

z! + v ’ + xz” - 


xz 12 df , , z' 


/3 o r ,/2 

z) + nrww +-= 0 


(4.20) 


Z 2f(z)dz y ' f(z) z 

In the presence of the B 2 field, the boundary condition on the probe D3-brane becomes 
the mixed boundary condition: 


(4.21) 

z=zo 

Here a is a dimensionless electric field defined in (j4.13j) . 

It is difficult to solve analytically the differential equation (I4.20p with the boundary 
condition (14.21j) . Hence let us solve it numerically and study the behavior of z(x). As a 
result, we will show the existence of two critical electric fields and the consistency of the 
values with the ones obtained from the potential analysis in subsec. 4.2. Hereafter, we 
take that 

A = 100 (2v tT f L 2 = 10), 



to validate the holographic description (i.e., large A). 

For some values of Zo/z t , numerical results are shown in Fig. [6] The left figure shows 
exponential suppression factors, and the right one is plots of the classical action. Typically, 
the suppression factors tend to vanish below certain values of a (i.e., E = E c ), depending 
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on values of z^/zt- This is the same ^-dependence as in the Coulomb phase. In the 
case of the confining phase, one can see a new behavior that the classical action diverges 
(i.e. the exponential factor vanishes) at the value E = E s . This result indicates that the 
Schwinger effect does not occur when E < E s . Note that the value of E s agrees with the 
potential analysis as well. As a matter of course, the value of E s depends on Zq/z^, and 
E s becomes zero as z t —> oo . 

In total, the two critical electric fields E s and E c have been seen from Fig. [ 6 ] These 
are the same as the ones obtained by the potential analysis in subsec. 4.2, 

L 2 


L 2 

E c = Tp—, 


E s = T F — 5 -. 


(4.22) 


This agreement supports that our numerical results of the production rate are consistent 
with the previous potential analysis. 




a 


(a) the exp factor 


(b) the classical action 


Figure 6 : The plots of the exponential factor and the classical action. 


The critical behaviors 

Let us examine the critical behaviors of the classical action numerically. We first argue the 
behavior of the classical action around E — E s , by focusing upon the singular behavior. 
Then the critical behavior around E — E c is discussed numerically and analytically. For 
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Figure 7: The behaviors of S Bs and C Bs (a s ) near E = E s . 


both critical behaviors, we will introduce new exponents and determine the numerical 
values. 


i) the critical behavior around E = E s 

Let us consider the behavior near E = E s . The log-log plot of Sb s shown in Fig. [7] (a) 
indicates that the a-dependence of Sb 2 may be described by 

= Cbs ( as ) ^ + the regular , a s = ^ . (4.23) 

(a - a s y E c 

Figure [7](b) shows that the coefficient Cb 2 (o : s ) is well approximated by 

CbM,) = ~( (4-24) 

As a result, by combining (j4.19j) and (I4.23p . the behavior of the Wilson loop radius R 
near E = E s is evaluated as 

R = -- ° + the regular. (4-25) 

a — a s 

This result means that R tends to diverge as E —* E s + 0. 

The next is to examine the NG action. The log-log plot shown in Fig. [8] (a) indicates 

that the a-dependence of the NG action is approximated by 

C NG (a s ) a T ) ng(«s) , , (A 

6 N g = -yr 4-4 the regular. (4.26) 

(a — a s ) a — a s 


27 
































(a) the leading term 


(b) the sub-leading term 


Figure 8: The leading and sub-leading terms of S'ng near E = E s . 


Here Cng(«s) and Dng(o!s) are scalar functions to be determined. 

We should comment on the validity of our numerical results. The plots in Fig. [S] would 
be valid for a s > 0.1 because of the limitation of numerical precision. If the accuracy of 
numerical analysis could be increased, then the lines in Fig. [8] are extended more to the 
right-hand side. 

In total, the total classical action S = S'ng + Sb 2 is given by 


5= c(a - )a „ + TM + theregu i ar , 


(a — a s ) 2 a — a, 
where we have defined new quantities as 


(4.27) 


C(a s ) = Cng(«s) + Cb 2 (&s) , D(a s ) = .Dng^s) 


Note that C(a s ) is nonzero for a s ^ 0. Therefore, in the deconhning limit a s —> 0, the 
total action S becomes 

rvni 4- mm 

(4.28) 


S= C (0 ) + O (°) +t he regular. 


a 


Then the following relation should be satisfied, 


C(0) + D(0) = ^, 


However, we have not confirmed this point yet because of the limitation of numerical 
analysis as mentioned in the last paragraph. 











It would be nice to consider the meaning of the power of the singularity. It is associ¬ 
ated with the critical behavior and hence may be interpreted as a new critical exponent. 
Suppose the following divergent form around E = E s : 

S = A(a s ){a - a s )~^ + ■ ■ ■ . (4.29) 

Now 7 S is a new exponent. Our numerical computation indicates that 

7s = 2. 

This may be regarded as a non-trivial prediction to the real QCD if this value could be 
interpolated to QCD via the universality argument. In this sense, it would be significant 
to check the universality of this exponent for general confining backgrounds. 

ii) the critical behavior around E = E c 

Let us next discuss the behavior near E = E c . The classical action should vanish as 
a —y 1, and hence one can expect the following behavior: 

S = B(a s )(l-ay c + ■■■ . (4.30) 

Here y c is a positive constant to be determined and B(a s ) is an unknown function. The 
log-log plot of the classical action is shown in Fig. [9] (a). The result indicates that the 
exponent y c is given by 

7c = 2. 

Moreover, from Fig. El (b), one can determine the asymptotic form of B(a s ) in the a s —> 0 
limit. As a result, the resulting classical action has been determined as 

S = Lhl - a) 2 + 0((1 - a) 3 ) . 

Note that this expression completely agrees with the asymptotic form derived from (13.201) 
analytically. This agreement ensures the consistency of our numerical computations. 

5 Conclusion and Discussion 

We have reviewed the recent progress on a holographic description of the Schwinger effect 
in the Coulomb phase and the confining phase. The part of the Coulomb phase is mainly 
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Figure 9: The behavior of the classical action near E = E c . The log-log plot in Fig. (a) 
indicates that 7 C = 2 universally. Figure (b) indicates that the coefficient approaches 5 
(= asa s -)-0 (Recall that A = 100). 

based on the seminal work by Semenoff and Zarembo [133. The part of the confining 
phase is a summery of a series of our works [T4l fTO] . 

First, we have introduced the world-line instanton method, which is a standard method 
to compute the Schwinger production rate. A remarkable point on this method is that the 
rate can be evaluated even at arbitrary coupling. Then a circular Wilson loop appears in 
the middle of the computation, in comparison to the weak-coupling analysis. This method 
is also applicable to the Af = 4 SYM theory through the Higgs mechanism by assuming 
the heavy quark mass and the weak-field condition. Notably, in this case, a circular 1/2 
BPS Wilson loop appears in the formula of the production rate (13.41) . The VEV of the 
Wilson loop is estimated by using the AdS/CFT correspondence. However, there are 
unsatisfactory points: 1) it seems unlikely that the Schwinger effect occur for very heavy 
quarks, 2) disagreement of the critical electric fields computed from the production rate 
and the DBI action of a probe D3-brane. 

To resolve this problem, Semenoff and Zarembo improved the AdS/CFT set-up. The 
probe D3-brane is put at an intermediate position in the bulk AdSs. Then the quark 
mass becomes arbitrary and the critical electric fields nicely agree. According to this 
improvement, the production rate is corrected with an additional term, which can be 
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ignored under the weak-field condition. Note also that the Coulomb potential is modified 
m as well. The critical electric field obtained from the potential analysis also agrees 
with the results from the production rate and the DBI action. 

The holographic method argued in the Coulomb phase has been generalized to con¬ 
fining gauge theories. In this review we have focused upon an AdS 5 soliton background. 
The potential analysis has been done [2] and it has been shown that there are two kinds 
of critical behaviors around (1) E = E s and (2) E = E c . When E < E s , the Schwinger 
effect cannot occur due to the confining string tension. In the region with E s < E < E c , 
it is possible as a tunneling effect as usual. When E > E c , the vacuum is unstable catas¬ 
trophically. Finally, the production rate has been evaluated numerically. Then we have 
introduced new critical exponents associated with the two critical behaviors. The values 
of the exponents are non-trivial results and might be regarded as a prediction in the real 
QCD via the universality argument, as in the case of the ratio of the shear viscosity r] to 
the entropy density s, r\fs = 1/47T [55]. Thus it is of importance to check the universal¬ 
ity of them for various backgrounds mt As a matter of course, it should be significant 
to reveal the mathematical foundation for the universality by employing the holographic 
approach discussed here as a compass. 

We have concentrated on an AdS soliton background as a confining geometry. As 
other examples, one may consider AT = 2 supersymmetric QCD and the Sakai-Sugimoto 
model [H|. The Schwinger effect in the former case has been studied in a series of 
works [471 - 149] . As for the latter case, see the recent work |50]• Another holographic 
Schwinger effect based on the bottom up approach is discussed by Dietrich [5Tj . It would 
be significant to study the Schwinger effect in confining theories from various perspectives. 
It is also interesting to study a holographic description of the Schwinger effect in de Sitter 
space [52]. It is argued that no pair production occurs in a plane-wave background [53] . 

We believe that the Schwinger effect in confining gauge theories would be a key in¬ 
gredient in looking for new aspects of QCD in the presence of extremely strong external 
fields. 
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